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"Radiative Transfer Models for the Earth Radiation Budget Studies." The 
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NAG-1 -21. The grant was monitored by Dr. John T. Suttles of the Radia- 
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ACCURATE EVALUATION OF HOMOGENEOUS 
AND NONHOMOGENEOUS GAS EMISSIVITIES 

By 

S. N. Tiwari 1 and K. P. Lee 2 
SUMMARY 

Spectral transmittance and total band absorptance of selected in- 
frared bands of carbon dioxide and water vapor are calculated by using 
the line-by-line and quasi -random band models and these are compared 
with available experimental results to establish the validity of the 
quasi-random band model. Various wide -band model correlations (such as 
Edwards and Balakrishnan, Tien and Lowder, and Felske and Tien) are 
employed to calculate the total band absorptance and total emissivity of 
these two gases under homogeneous and nonhomogeneous conditions. These 
results are compared with available experimental results under identical 
conditions. From these comparisons, it is found that the quasi-random 
band model can provide quite accurate results and is quite suitable for 
most atmospheric applications. Edwards and Balakrishnan, and Tien and 

Lowder correlations can provide reasonable results for homogeneous and 
♦ 

nonhomogeneous conditions and, therefore, are useful for most engineer- 
ing applications. 

^•Eminent Professor, Department of Mechanical Engineering and Mechanics, 
Old Dominion University, Norfolk, Virginia 23508. 

Graduate Research Assistant, Department of Mechanical, Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 



Chapter 1 


INTRODUCTION 

The concept of total emissivity has proven to be very important in 
many radiative transfer analyses. It is useful in calculating atmos- 
pheric radiation fluxes and cooling rates [1-3]*, and predicting the 

emission and absorption of infrared radiation by hot gases in the heat 

» 

transfer within high temperature furnaces, combustion chambers, heat 
engines, and chemical reactors [4-8]. Moreover, it is very convenient 
to use the emissivity charts of Hottel and co-workers [9] to estimate 
heat transfer from homogeneous, as well as nonhomogeneous, gases quite 
accurately. But, in tile extrapolated portions of the emissivity charts, 
especially for water vapor, a discrepancy exists [10, 11]. In order to 
examine the usefulness of Hottel 's data, a number of investigators cal- 
culate the emissiviti'es which are based on spectral data or acquire them 
from experiments to make a comparison [4, 5, 10, 12-17]. 

The most accurate theoretical procedure to compute the spectral gas 
properties of a vibration-rotation band is probably the direct integra- 
tion (line-by-line) method, which calculates the monochromatic absorp- 
tion spectrum by numerically summing the contributions of the individual 
molecular absorption lines. Although this method is of high spectral 


*The numbers in the brackets indicate references. 
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resolution, it is time-consuming and requires knowledge of ..idividual 
line parameters. If one wants to save computational time and does not 
require high resolution, one can use narrow band models to evaluate the 
spectral absorptance within a narrow frequency interval of a vibration- 
rotation band. The commonly used narrow band models are El sasser, sta- 
tistical, random-El sasser, and quasi-random models. The quasi-random 
band model is probably the best one to represent emission of a vibra- 
tion-rotation band accurately. 

Although the computational time required when using the narrow band 
model is much less than the time required when using the line-by-line 
method, it still requires a large amount of computer time, as well as a 
massive amount of input data, which makes it impractical for engineering 
applications. Therefore, it is desirable to make use of the so-called 
"wide band models" [11, 15, 18-21] which represent absorption from an 
actual band with reasonable accuracy. As these models are based on 
various narrow band model relations for absorption by assuming that the 
line intensity is an exponentially decaying function of the wavenumber, 
they are called the exponential wide band models. In addition to reduc- 
ing the mathematical complexities in radiative transport calculations 
and saving the computational time, it is convenient to express the total 
band absorptance or the band emissivity by fairly accurate continuous 
correlations [22-25]. 

Assuming the real optical paths of combustion chambers, furnaces, 
or the atmosphere to be homogeneous can lead to considerable errors when 
an estimate is made of the total radiative flux. Since optical paths 
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usually are nonhomogeneous, actual radiative flux calculations are very 
complicated and time-consuming. In order to save a considerable amount 
of computational time and make the problems more simple and suitable for 
engineering applications, one can apply the Curti s-Godson scaling ap- 
proximation to predict the band emission from the nonhomogeneous gases. 

The purpose of this study is to calculate the spectral and total 
emissivity of carbon dioxide and water vapor by using band models and 
band model correlations. We specify these two gases because they are 
the most effective in the emission and absorption of infrared radiation. 
The spectral transmittance and total band absorptance of different bands 
of these two gases, as obtained by using the quasi -random band model and 
several band model correlations, are compared with available experiment- 
al results [26-28] in order to establish the validity of a particular 
formulation to specific applications. 

The basic theoretical formulations of emissivity for homogeneous 
and nonhomogeneous gases are presented in Chap. 2. A brief review of 
line and band models, and band model correlations are given in Chap. 3. 
Physical conditions and data sources for this study are presented in 
Chap. 4. Results are presented and discussed in Chap. 5. 


Chapter 2 


THEORETICAL FORMULATION 

A large amount of information concerning various relations for 
emissivity of different gases is available in the literature [1-2, 10- 
15, 29-31]. In t’.iis section, we only review the basic formulations for 
calculating the emissivity for homogeneous and nonhomogeneous paths 
briefly. 


2.1 Homogeneous Gas Emissivities 
For a homogeneous path of absorbing-emitting gases, the spectral 
emissivity, e^, at wavenumber w is defined as 

e a 3 A u 3 i-V = 1-expf-K^X) (2.1) 

where A w is the spectral absorptance, is the spectral transmit- 
tance, is the volumetric absorption coefficient and X=py is the 
pressure path length. The total emissivity, e^, is obtained by inte- 
grating over the entire energy spectrum; it is defined by 

00 

e T (X,P,T) = — [ e M (X,P,T)B M (T)dw (2.2) 
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where a is the Stefan-Boltzmann constant and B, is the Planck func- 
tion evaluated at wavenumber to. if one considers only the band emis- 
sivity, of the nth band, Eq. (2.2) is integrated over the band 

region only. 

For a single gas, where the absorption spectrun consists of vibra- 
tional-rotational bands, the total emissivity may be expressed as 

e t (X,P,T 2 B.(T)A.(X,P,T) (2.3) 

1 aT 1 * i 1 1 

where B^(T) is the Planck function evaluated at the band center, and 
A. is the integrated (total) band absorptance of the ith band and is 
given by 


A = / A^do) = / [l-exp(-K aj X)]dw (2.4) 

Ao) tun 

By using an appropriate band absorptance relation, the total emis- 
sivity can Be calculated with Eq. (2.2) and Eq. (2.3). Especially, one 
can obtain the emissivity without the integration by employing the band 
model correlations; these are extremely convenient and useful for engi- 
neering applications. 

The spectral emissivity for a two component mixture of radiating 
gases denoted by subscripts 1 and 2 is given by 

e w,12 s Vl V2 s l-exp[-(K (|1> i x 1 + K Ui2 X 2 )] 


(2.5) 


Thus, the total emissivity is expressed as 


e T,i2 ■ t e M2 (X>P - T)B <»< T)d “ 

or 

g T,12 ■ l B «. (T )1 1 -«P[-< K 0J ,1 X 1 + K <.,2 X 2 ): » d “ < 2 - 6) 

Equation (2.6) may be rewritten in the following form: 

e T,12 = e T,l + e T,2~ 

00 

JL. [ B^tTJCl-expt-K^ iXjJJCl-expC-K^ 2 x 2 )3<jm (2.7) 

where Cy ^ and ^ are the total emlssivities of the pure compo- 
nents 1 and* 2, respectively. The last term of Eq. (2.7) is called the 
correction for partial overlapping which is discussed in detail by 
Penner and Varanasi [32]. 

2.2 flonhomogeneous Gas Emissivities 
For a nonhomogeneous path, temperature, pressure and concentration 
of gas vary along the path and this affects the absorption of the gas; 
and, therefore, the emissivity relation should be modified appropriate- 
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ly. Corresponding to Eq. (2.2), the emissivity for a nonhomogeneous 
path for the path between level x and x* can be expressed by 

<*» B, fx') 

e(x,x')=/ A (x,x') -- — dw (2.8) 

' w B(x') 

where A 6J (x,x') is spectral absorptance between x and x', and 
B(x') = /B aJ (x')dw = crfV*. Thus, it is necessary to acquire the 
spectral absorption A^x.x') before evaluating the emissivity. 

Nonhomogeneous behavior can be considered by employing the Curtis- 
Godson approximation [1], which appropriately averages the mean line 
intensity and line spacing over the nonhomogeneous path. For noniso- 
thermal gases, this approximation, by using narrow band models, is not 
suitable for most engineering applications because of the complexity of 
the models and the large amount of computational time required. It is, 
therefore, desirable to use wide band models. Several investigators 
[19, 33-36] have applied the Curtis-Godson approximation to wide band 
models by obtaining three scaled parameters for the entire band. These 
are: a scaled integrated intensity, a scaled line width to spacing 

ratio, and a scaled band width parameter related to the rotational con- 
stant. 

The simplest and most accurate model is the one suggested by Chan 
and Tien [33]. For this model, the scaling parameters can be expressed 
as 
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L 


(«) h = (1/x) 

/ apdy 

(2.9) 


L 


(eg) h = (l/« h x) 

/ 3P apdy 

(2.10) 

(u) h = (l/a h x) 

l 

/ awpdy 

(2.11) 


L 

where x » / p dy, L is the thickness of the gas layer, p is the 
© 

density, y is the distance along the optical path; a, u> and 3 are 
integrated intensity, line width to spacing ratio parameter and band 
width parameter respectively; and P fi is the effective broadening 
pressure. By using these three scaled parameters in an available band 
model correlation, one can determine the nonisothermal band absorptance 
and band emission which is given by 

x l_ . 

E = / B CT (x) ] — A[T(x' ),x]dx (2.12) 

• dx 

By employing Eq. (2.12), one may obtain the nonhomogeneous gas emis- 
sivity. 


Chapter 3 


BAND MODELS AND CORRELATIONS 

Various theoretical formulations of spectral, narrow and wide band 
models, as well as band models correlations, are available in the liter- 
ature [1, 11, 18, 20-25, 37-43]. In this section, these models and cor- 
relations are described briefly. 

3.1 Direct Integration (Line-By-Line) Model 

In this method, the entire spectral range of interest is first 
divided into a large nunber of narrow intervals. Each interval is then 
divided into a variable number of subintervals depending upon the number 
of lines within the interval. Two very narrow subintervals are created 
on each side of a line center. The transmittance, and then the spectral 
emissivity, is computed at four frequency locations in each subinterval 
and Is finally averaged over each interval. 

The total absorption coefficient in Eq. (2.1) at wavenumber to is 
computed into two parts as 

K(to) * K D (<o) + K^(to) s (3.1) 

/ 

D W 

where K (to) and K (to) are called the direct and wing contributions 
respectively. Direct contribution originates from lines in very close 
vicinity (on both sides) and is obtained from 
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wt. 


K°(w) = £ s n Y n /^ 7r [( £0 " 0> n ) 2 * Y n^ (3.2) 

where oi^ refers to the nth contributing line, S n and Y n are the 
Intensity and the Lorentz line half-width respectively. The wing con- 
tribution arises from lines which are further from to than the range of 
direct contribution and is obtained from 

K W (oj) = £ S n Y n /[Tr(<o-to n ) 2 ] (3.3) 

Both Eq. (3.2) and Eq. (3.3) are based on the assumption that the line 
shape is Lorentzian. For a complete discussion on the direct integra- 
tion procedure, one should refer to [38-40], 

3.2 Narrow Band Models 

The mean spectral emissivity within a narrow frequency interval can 
be represented with reasonable accuracy by the so-called "narrow band 
model." Four commonly used narrow band models are Elsasser, statisti- 
cal, random Elsasser, and quasi-random. The application of any model to 
a particular case depends upon the nature of the absorbing-emitting 
molecule. Complete information on these models are available in the 
literature [X, 11, 24]. 

3.2.1 Elsasser (Regular) Band Model 

The regular Elsasser band model consists of equally spaced Lorentz 
lines of equal half-width and intensity [44]. The absorption coeffi- 
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clent Is a periodic function (with the period of the line spacing) and 
is given by 


( 1 ) 


oo 

= E 
n=-» 



t j 

(w-nd) 2 + Yj 


(3.4) 


where w is the distance from the center of any line and d is the 
distance between adjacent lines. Elsasser [44] showed that Eq. (3.4) 
can be expressed in an alternate form as 


where 


= !L sinhp 
d coshu - cosZ 


(3.5) 


e = 2iTYj/d, 


and Z = 2ir(w-a» o )/d 


(3.6) 


The mean spectral emissivity of the periodic line pattern over the line 
spacing d is obtained by combining Eqs. (3.5) and (2.1) as 

e = 1 - -1 /* exp[ ]dZ (3.7) 

2ir “ir cosh0 - cosZ 

where 


x => SjX/2tt Yj 


(3.8) 


The approximate forms of the Elsasser model can be expressed as: 

1. The weak-line approximation 

e a l-exp(-ex) 0 + «, x«l (3.9) 

2. The strong-line approximation 

e =erf[{i0 2 x) ] x»l (3.10) 

2 

where erf(t) denotes the error function. • 

3. The nonoverlapping-line approximation 

e = 0L(x) 0«1, 0 2 x«l (3.11) 

where L(x) Is Ladenberg-Reiche function defined as 
♦ 

L(x) = xexp(-x) (I Q (x) + I x (x)] (3.12) 

where I Q and 1^ are the Bessel functions of imaginary arguments. 

4. The strong nonoverlapping-line approximation 

, 1/2 

e = 0 (2 x/tt) x »1, 0«1, 0 2 x«l (3.13) 


3.2.2 Statistical (Mayer-Goody) Model 

The statistical band model is based upon the assumption that, in a 


given wavenunber interval, the spectral lines are spaced randomly, and 

the intensity of these lines can be specified by some distribution func- 

> 

tion [1]. 

For uniform statistical model, which assumes all spectral lines in 
a narrow spectral interval to be equally Intense, the mean spectral 
emissivity for lines of the Lorentz profile is expressed as 

e = l-exp[-0xe" x (l o <x) + ^(x))] (3.14) 

The line spacing d in the expression p = (2tr Yj/d) refers to average 
spacing over the narrow band. The approximate forms of this model are 
the same as that of the El sasser model except the strong-line approxi- 
mation which is expressed as 

1/2 

e = l-exp[-(20 2 x/ir) ] x»l (3.15) 

These should be expected because the particular arrangement of the spec- 
tral lines in the band does not influence the emission in the weak-line 
approximation and only intensity distribution function influences the 
emission in the nonoverlapping-line approximation. 

The general statistical model assienes the distribution of line in- 
tensities to be, generally, represented by an exponential distribution 
(i.e., in a band, the probability of finding a spectral line of inten- 
sity Sj in a given intensity range decreases exponentially). Then the 
mean spectral emissivity by this model consisting of Lorentz lines is 



I 


e = l-exp[-px/(l + 2x) ] 


(3.1G) 


Now the line spacing d and line Intensity S in $ and x should be 
replaced with the average value over the narrow spectral Interval. The 
weak-line and strong-line approximation are the same as that in the uni- 
form statistical model, but the nonoverlapping-line approximation for 
this model is expressed as 


e + 3x/(l + 2x) 


(3«1, e 2 x«l 


(3.17) 


Other line intensity distributions also have been suggested for the 
statistical band formulation; for detailed discussion on this, one 
should refer to [24, 25]. 


3.2.3 Random El sasser Band Model 


The random El sasser band moael assumes the random superposition of 
several different El sasser bands. Each of the superposed bands may have 
different line intensities, half-widths and spacing. As many different 
El sasser bands as necessary may be superimposed in this model. For N 
randomly superposed Elsasser bands, the mean spectral emissivity is 
obtained by the relation 


e-l-n [1-A r (x<, 


(3.18) 
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/v 

where A^ (x^, 3.j) is the absorptance of the 1th Elsasser band and is 
given by 


A E,1 /S i " l A E,1< x 1- 0 1 )P E< S i' S o1> dS i < 3 ' 19 > 

where P^(S^ ■ S o1> 15 the line intensity distribution function. 

If an exponential distribution of intensity for narrow Elsasser 
bands is assumed, then 

P E (S i’ S oi ) * (1/S oi> exp(-S i /S oi ) (3.20) 

and Eq. (3.19) becomes 

A E,/ 6 i = ^i x oi s1n h0 i )/C(3 i x oi sinh6 i + ccshB .) 2 -l] 1/2 (3.21) 

♦ 

A combination of Eqs. (3.18) and (3,21) yields the expression for mean 
spectral emissivity by a modified random Elsasser band model as 

i/2 

e = (3x Q sinh3)/[(3x 0 sinh3 + cosh3) 2 -l] (3.22) 

The weak-line and strong-line approximations of the random Elsasser 
model can be expressed, respectively, as 


B + X«1 


(3.23) 


e = 1- n [exp -B-xJ] 

1=1 1 1 

and 

N 1 1/2 

e = 1- H {l-erf[( i B^xJ ]} x»l (3.24) 

1=1 2 1 1 

The approximate forms of the modified random El sasser model are given by 
the same relations as the general statistical model except the nonover- 
lapping-lina approximation which depends on the line intensity distri- 
bution assuned. 


3.2.4 Quasi -random Band Model 

Quasi-random band model assumes neither a regular nor a random 
spacing of the spectral lines. The entire spectrum is divided into a 
number of narrow frequency intervals with equal spectral width. The 
lines within each narrow interval are as c, Jiied to be distributed randomly 
and grouped into five intensity decades. Average transmittance for the 
interval is computed for lines in each decade separately and final 
transmittance is obtained by multiplying them for the five decades. 

Then, using Eq. (2.1), the mean spectral emissivity is obtained. 

The average transmittance over interval, 6, due to a single line 
may be expressed as 


T n (o») = - / exp[-S n Xf(w,w n )]dti) n 
6 6 


(3.25) 
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where f(w, to^) is the shape factor for the line with center at w 0 . 

The average transmittance, due to all the lines in that decade is 

given by 

i - N 

T d M = { - / exp[-S n Xf(o),<o n )]d(o n } (3.26) 

where N is the number of lines and $ n is the average intensity of 
all the lines within the decade. At last, the average transmittance due 
to all lines in the five intensity decades of the interval is expressed 
as 


5 5 ~ N 

T k ( '“> * ", T d ( “) = ", i — / exp[-S n Xf ((.,«> )]du } (3.27) 

d-1 d=l 6^ 5^ 

where k represents the kth interval of entire spectrum. The wings of 
the lines in the adjacent intervals also make a significant contribution 
to the absorption in 5^. Therefore, the resultant transmittance over 
the interval 6^ is given by 


k 

T k<“> * Vk<“> " T k-j(“) (3-28) 

j*k 

T k-j{ w } re P rese nts the transmittance in <5^ due to lines in 


where 
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5 j . For a complete discussion on the quasi-random band model, one 
should refer to [37, 39, 41]. 

The model formulations discussed above represent the emissions of 
gases under homogeneous conditions. They can be extended to nonisother- 
mal nonhomogeneous optical paths by employing various scaling approxi- 
mations available in the literatures [1, 19, 33-36, 39, 46, 47]. 

The assumption of either regular or random spacing for spectral 
lines may be a serious defect of the Elsasser and statistical models 
because the spacing between the spectral line is in general neither 
■uniform nor random. Furthermore, the wing contribution cannot accurate- 
ly be taken into account by these models. Therefore, the quasi -random 
band model is probably the best model to represent the emission of a 
vibration-rotation band. In this study, only the quasi-random band 
model is used to make comparison with other formulations. 

3.3 Wide Band Models. 

Wide band models can be used to predict band emissivity without a 
lot of calculations but with reasonable accuracy. They are very simple 
and convenient for engineering applications. Detailed discussions on 
the wide band models are given in [11, 15, 18, 20, 24], 

3.3.1 The Box Model 

This model, first introduced by Penner [12], assumes that the 
absorption coefficient, K^, is constant over an effective band width 
Aw. The expression for the total absorptance by this model is given by 
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A = / [l-exp(-K w X)]dw = (A“>) e [l-exp(-KX)] (3.29) 

where (Aid) is the effective band width and K is the mean absorption 
6 

coefficient for the interval (Aw) e . 

3.3.2 The Exponential Wide Band Model 

Edwards et al. [18, 19] have considered various wide band models 
for the absorption of vibrational-rotational bands and have concluded 
that three parameters, which are the mean line intensity to spacing 
ratio. S/d, the mean line width to spacing ratio, iry/d, and the 
effective broadening pressure, P Q , are necessary for a complete de- 
scription of band absorption. The final forms of the total band 
absorptance relation presented by Edwards and Menard [19] and Edwards 
and Balakrishnan [21] are based on the formulation of the narrow sta- 
tistical band model; and these are stmmarized in Table 3.1 and 3.2. 

The other formulations of total band absorption with different 
narrow band models are discussed in detail in [24]. 

3.4 Band Model Correlations 

In order to reduce the mathematical complexities and save the com- 
putational time, it is often desirable to express the integral form of 
the total band absorptance or the band emissivity by fairly accurate 
continuous correlations. Several continuous correlations for the total 
absorptance of a wide band are available in the literature. These are 
briefly discussed here in the sequence that they became available in the 
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Cl, C 2 , C 3 , b and n are correlation parameters [11,19] , X is mass patb length, g/m 2 , 
Pe-t^k + *o = ^ a*®, absorbing-gas partial pressure, is the 
broadening-gas partial pressure. 
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literature. 

3.4.1 Tien and Lowder Correlation 

In 1966, a continuous correlation was proposed by Tien and Lowder 
[22], and this is of the form 

A(u»p) = A{u,0)/w 

= ln(uf(t){(u+2)/[u+2f(t)]}+l) 

where 

f (t) = 2.94[l-exp(-2.6t)], t=P/2 

A(u,&) is dimensionless total band absorptance, u=SX/w is 
dimensional path length and 0 = 2iry/d is the line structure 
This correlation does not reduce to the correct limiting form 

rv 

square-root limit [20], and its use should be made for p>0.1. 

3.4.2 Goody and Belton Correlation 

In 1967, another continuous correlation was proposed by Goody and 
Belton as [48] 


{3.30) 


(3.31) 

the non- 
parameter, 
in the 


" 1/2 

A(u,e)=21n{l+u/[4+jru/4t)] } 


(3.32) 
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Although this correlation satisfies the linear, square-root, and loga- 
rithmic limits, its use is restricted to relatively small 0 values. 

3.4.3 Tien and Ling Correlation 

In 1969, a simple two-parameter correlation was developed by Tien 
and Ling as [49] 


A(u) = sinh X {u) (3.33) 

This is valid only for the limit of large 0. 

3.4.4 Cess and Tiwari Correlations 

A relatively simple correlation has been proposed by Cess and 
Tiwari as [20] 


” " 1/2 
A(u,0) = 21 n (1+u/ { 2+[u (1+1/0 ) ] }) 


(3.34) 


where 0 = 4t/ir. The use of this correlation is justified at relatively 
high pressure for gases whose spectral behavior can be described by the 
general statistical model. By si ightly modifying Eq. (3.34), another 
form of the wide band absorptance was obtained as [24] 


" 1/2 
A( u ,0 ) =21 n ( l+u/{ 2+[ u( C+n /20 ) ] }) 


(3.35) 


0.1, 6< and ail u values 

0.1, 0> 1 and i* < 1 

0.25, 3> 1 and u > 1 


Both Eq. (3.34) and Eq. (3.35) reduce to all the limiting forms [20, 


3.4.5 Edwards and Balakrishnan Correlation 


In 1972, a continuous correlation was proposed by Edwards and 
Balakrishnan as [50] 


A(u) = 1 n( u) + E x (u) + Y + — - E 3 (u) 

2 


This is valid for large pressure and at large path lengths. 


3.4.6 Felske and Tien Correlation 


In 1974, Felske and Tien proposed a correlation for A(u,B) as 


A(u,0) = 2Ej, (tp u ) + Ei(p u /2) - Ei[(P u /2)(l+2t)] + 


ln[(tp..) 2 /(l+2t)] + 2y 


[Q 





p u = { (t/u) [l+(t/u)]J 


(3.38) 


Ei ( x) is exponential integral of the first order, and Y=0. 5772156* is 
the Euler's constant. This correlation is based on the general statis- 
tical model and its validity is claimed for the entire range of the 
governing parameters. 


3.4.7 Tiwari and Batki Correlation 

In 1975, Tiwari and Batki proposed the following correlation [24, 


A(u) = y + ln(u) + Ei ( u) 


(3.39) 


This correlation is valid for all path lengths but for t=(0/2)>l. 


3.4.8 Wang-Correlatiori 

A continuous correlation has been proposed recently by Wang as [25] 


A(z) = exp(4^i ) {Ei [20! (1+z) ] - (43i) } + In (1+z)- 


ln 2 + Ei [20i (z-1)] + In [2Bi(z-l)] + y 


(3.40) 


where Bi = e/4, z=(l+u/Pi) . Thi s total band absorptance correla- 
tion is obtained by employing the Malkmus narrow band model [45] and the 


■ 

B 
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Edwards and Menard wide band approximation [16]. The predicted total 
band absorptance of this correlation is always smaller than the Felske 
and Tien. 

Band absorptance results of various correlations are compared and 
discussed in some detail in [24, 51]. The Felske and Tien correlation 
agrees well with the general statistical model while the Tien and Lowder 
correlation is in general agreement with the El sasser model. In this 
study, only these two correlations are employed to calculate the total 
emissivity. 


Chapter 4 


PHYSICAL CONDITIONS AND DATA SOURCES 

The primary atmospheric absorbers and the primary radiating species 
of the hot gases of combustion are the carbon dioxide and water vapor 
due to their comparatively high absorptlvities and emissivities in the 
near Infrared region. Tbe total emissivities of these two gases were 
calculated by using the quasi -random band model {only at room temper- 
ature), Edwards and Balakrishnan wide band model relation, Tien and 
Lowder correlation and Felske and Tien correlation under conditions of 
temperatures at 300, 600, 1000, 2000 K, pressures at 0.1, 0.5, 1.0, 5.0, 
10.0 atm, and path lengths from 0.1 to 1000 atm-cm. Moreover, results 
of transmittance and total band absorptance were obtained for several 
bands of these two gases under conditions for which experimental meas- 
urements were available. 

The lipe parameters (position, strength, line width etc.) needed 
for this study were obtained from McClatchey c- al. [52]; the 
"McClatchey Tape" is available at the NASA Langley Research Center. Ro- 
tational and vibrational partition functions required to account for the 
temperature dependence of the line strengths were taken from McClatchey 
et al. [52]. Table 4.1 lists the significant bands of CO 2 and H 2 O 
which provide the essential contributions to the gas emissivity. 

The band absorptance parameters (the integrated band intensity, a, 
the line-width parameter, 0 , and the band width parameter, o») for the 
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Table 4.1 Significant band regions 
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wide band models and correlations are obtained from Edwards et al. [21, 

15]. Table 4.2 lists these parameters for CO 2 and H 2 O. The 
temperature and pressure dependence of these parameters are given by the 
following relations: 


o(T) = ct 0 


[l-expf-sJL^ u k fi k )MT) 
[l-expf-^ u 0>k s k )MT 0 ) 


(4.1) 


B(T) = no(T/T„) «(T)/ t (T 0 ) 


(4.2) 


“(T) = « 0 (T/T 0 ) 


(4.3) 


where 


?(T) = 


“ (v. + g.+ 5.-1)! -u.v t 

jj jj AAA q IS K. 

k< Wk <V 1)! v k‘ 

™ z <v K* g K~ 1)! e ” u k v k 
k=1 V 0 '3k- 1 ) ! v k ! 


(4.4) 


*(T) = 


m « 

n z 

1(31 Wk 

m » 

n z 

k=1 v k =v 0,k. 


e - u k v k l 1/z 

(9k- 1 ) 1 V = 

I’k^t * 8 ^ 1 e -Vk 

(9 k -D! v k ! 


(4.5) 


u k = he v k /kT, U 0>k = hcv k /kT c 


(4.6) 


in which T = 100 K, and hc/k £ 1.4388 cm-K 
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Chapter 5 

RESULTS AND DISCUSSIONS 

Spectral transmittance of important bands of carbon dioxide and 
water vapor are evaluated by employing the line-by-line (LBL) and quasi- 
random band (QRB) models; and these are compared with available experi- 
mental results. The total band absorptance and the total emissivity of 
these two gases are calculated by the QRB model, Edwards and 
Balakrishnan wide-band absorptance relations, Tien and Lowder correla- 
tion, and Felske and Tien correlation. The nonisothermal band absorp- 
tance and emission, however, are calculated only by using the wide-band 
model correlations. All these results are compared with experimental 
results available in the literature. 

5.1 Homogeneous Results for Carbor. Dioxide 

The spectral transmittance and total band absorptance results of 
selected bands (15y CO 2 , lty CO^, 4.3p C0 2 and 2.7y C0 2 ), as 
well as the results of total emissivity, of carbon dioxide are presented 
in this section. 

5.1.1 15m CO 2 Band 

The spectral transmittance results for this band, as calculated by 
the QRB model, are shown in Fig. 5.1 along with experimental results. 

The agreement between these results is quite good over the entire band 
pass. 
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Pig, 5.1 Comparison of transmitt ances of the 15}i COg band 
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The total band absorptance results were calculated by the QRB model 
{only at 300 K) and wlde-band model correlations at four different temp- 
eratures and are presented In Tables 5. 1-5. 4. At 300 K, the QRB model 
results agree very well with experimental results at large path lengths 
but give lower values than the experimental values for very small path 
lengths. Among the results of the correlations, Edwards and 
Balakrishnan, and Tien and Lowder correlations generally show better 
agreement for medium and high pressures as well as for large path 
lengths, while Felske and Tien's correlation yields much lower values 
than experimental values. For low pressures and small path lengths, 
however, Felske and Tien's correlation provides better agreement; 

Edwards and Balakrishnan, and Tien and Lowder correlations give higher 
values than experimental values. 

5.1.2 10|i CC^ (Hot) Band 


The spectral transmittance results, as calculated by the QRB and/or 
line-by-line models, are presented in Figs. 5.2 and 5.3 for two differ- 
ent pressures and path lengths. Hie QRB model results show very good 
agreement with the experimental results in Fig. 5.2; however, the QRB 
model results exhibit a slightly lower absorption in Fig. 5.3 for higher 
pressure and longer path length. This may be because the QRB formula- 
tion neglects the contributions of weaker lines under these conditions. 

The agreement between experimental and line-by-line results is seen to 
be excellent. 
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Table 5»3 Comparison of total band absorptance for the 
15p CO2 band at T=833 K 



PD=( (A-A e )/Ae)X100 


Table 5.4 Comparison of total band absorptance for the 
15>i COg band at T=1110 K 
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PD=( ( a-a e ) /Ae )X1 00 


1.0 
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5.1.3 4.3m COj Band 

The spectral transmittance results of the QRB model are compared 
with the experimental results in Fig. 5.4 for P 0 = 137 mm Hg and X = 
0.195 atm-cm. The agreement between the results is seen to be good 
except in the central portion of the band where the QRB model results 
exhibit a slightly lower absorption. 

The total band absorptance results for this band, as calculated by 
the QRB model (only at 300 K) and different correlations are compared 
with the experimental results in Tables 5. 5-5. 8 for four different temp- 
eratures. The QRB model results are seen to compare very well with the 
experimental results except for low pressures and small path lengths. 
Among the results of correlations, the Edwards and Balakrlshnan, and 
Tien and Lowder results show, in general, reasonable agreement with the 
experimental results; but the comparisons are not very good for low 
pressures and small path lengths. The Felske and Tien's results are 
lower for a]most all cases considered. 

5.1.4 2.7m CG^ Band 

The transmittance results of the QRB model are compared with the 
experimental results in Figs. 5.5 and 5.6 for different pressures and 
path lengths. In general, the agreement between these results is seen 
to be good, except that the QRB results show a consistently lower ab- 
sorption (higher transmittance) in Fig. 5.5. 

















PD=( ( A-A e ) /A b )X1 00 




-Pig. 5.5 Comparison of transraittances of the 2,7ji CO 2 band 





1.0 


ORIGINAL PAGE f3 
OF POOR QUALITY 



CO VO Tt* 

• • • 

o o o 


CM 

o 


o 

o 



VO 

in 


60 

•H 

ft 



aouviiirasimii 


48 


The total band absorptance results, as calculated by the QRB model 
(only at 300 K) and wide-band correlations, are presented in Tables 5.9- 
5.12 for four different temperatures. It should be noted that the QRB 
model results are not In good agreement with the experimental results 
for low pressures and small path lengths; however, they are in good 
agreement for large path lengths. Inspecting the results of wide-band 
correlations, it is noted that the Edwards and Balakrishnan, and Tien 
and Lowder results compare well with the experimental results except for 
a few cases considered. The Felske and Tien's correlation, however, 
always provides lower values. 

5.1.5 Total Emissivity for C0 2 


By using the QRB model (only at 300 K) and wide-band correlations , 
the total emissivity has been calculated according to Eq. (2.2) and Eq. 
(2.3) for the path lengths from 0.1 to 1000 atm-cm in order to make a 
comparison with the measured emissivities of HotteT et al. [9]. These 
results are presented in Figs. 5.7, 5.9, 5.11, and 5.13 for -four differ- 
ent temperatures and p * 1 atm. The QRB model results are seen to be 
in quite good agreement for all the path lengths. Among the correla- 
tions, the Edwards and Balakrishnan, and Tien and Lowder results show 
good agreement with the experimental results for long path lengths. The 
Felske and Tien’s correlation yields lower values for all path lengths 
in all cases considered. 

The total emissivity results obtained by using the wide-band corre- 
lations are compared in Figs. 5.8, 5.10, 5.12 and 5.14 for four differ- 


Table 5*9 Comparison of total band absorptance for the 
2.7}i CO 2 ttand at T=300 K 
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Table 5,10 Comparison of total band absorptance for the 
2.7}i C0 2 band at T=555 K 



Table 5.11 Comparison of total band absorptance for the 
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ent temperatures and five different effective broadening pressures. It 
is noted that the Felske and Tien's results are always lower than the 
results of other two correlations. This can be expected because the 
general statistical model always predicts lower absorption than the 
El sasser model due to more overlapping of the spectral lines in the 
statistical models. At low pressures, the Tien and Lowder's results do 
not follow the trend of other results. This would be expected because 
the low-pressure situation corresponds to the case of the square-root 
limit which is not satisfied by the Tien and Lowder's correlation. 
Therefore, use of the Tien and Lowder's correlation is not recommended 
for low pressures. For large path lengths, the Edwards and 
Balakrishnan, and Tien and Lowder correlations will give close results. 
For small path lengths, the Edwards and Balakrishnan' s correlation pro- 
vides lower results at low pressures while it yields higher results at 
high pressures than the Tien and Lowder's correlation. 

. 5.2 Homogeneous Results for Water Vapor 

The spectral transmittance and total band absorptance results of 
selected bands (rotational, 6.3u H 2 O, 2.7 y ^0 and 1.87 y H 2 0), as 
well as the results of total emissivity, of water vapor are presented in 
this section. 

5.2.1 Rotational Band 


The results for spectral emissivities obtained by the QRB model are 
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compared with the experimental results of Ludwig et al. [17] In Figs. 
5.15-5.20 for six different temperatures and path lengths. The QRB 
model results, In general, are seen to be higher than the experimental 
data. The results for T = 590 K (Fig. 5.15) appear to be relatively 
in good agreement Indicating that at low temperatures the QRB model will 
provide accurate results. The agreements, however, are not good for the 
temperature range of 800-1,800 K. For T 3 2,000 K, the QRB results 
are lower at shorter wavelengths and relatively higher at longer wave- 
lengths. The behavior exhibited in these figures could be due to a com- 
bination of several factors; but, overlapping of the spectral lines at 
higher temperatures will play the major role. At higher temperatures, 
the spectral lines overlap considerably and the QRB model does not 
account for this effect accurately. However, one can expect the QRB 
model to yield accurate results at temperatures closer to the room temp- 
erature since the line parameters are compiled at this temperature. 

5.2.2 6.3ii H 2 0 Band 


The comparison of spectral transmittance of this band is shown in 
Figs. 5.21a-5.21c. In order to present these results clearly, it is 
necessary to show them in three figures. The QRB model results show a 
quite good agreement with the experimental results. 

The total band absorptance results of this band are shown in Tables 
5.13-5.16 for four different temperatures. The QRB model results pre- 
sented in Table 5.13 for T = 300 K provide higher values at low 
pressures; however, the results agree very well with the experimental 
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results at high pressures and long path lengths. Anong the results of 
the correlations, it is noted that all the correlations underestimate 
the absorptance at low temperatures. At high temperatures, the Edwards 
and Balakrishnan, and Tien and Lowder correlations can predict the total 
band absorptance with reasonable accuracy. The Felske and Tien's re- 
sults are always found to be lower than the experimental results by 
about 20 to 40 percent. 

5.2.3 2.7m HgO Band 


Figures 5.22 and 5.23 illustrate a comparison of the spectral 
transmittance results, as calculated by the QRB model, with the experi- 
mental results for two pressures and path lengths. These results dis- 
agree in the wing regions of the band, although they agree well in the 
central portion. 

The total band absorptance results calculated by the QRB model and 
wide-band correlations are compared with the experimental results in 
Tables 5.17-5.20. The QRB model is found to estimate the absorptance of 
this band quite accurately. Anong the results of correlations, the 
Edwards and Balakrishnan, and Tien and Lowder correlations yield good 
agreement with experimental results for high pressures and large path 
lengths; however, they overestimate the band absorption by about 20 per- 
cent for low pressures and small path lengths. Note again, that the 
Felske and Tien's correlation usually underestimates the total band 
absorptance by about 20 percent. 











lable 5.17 Comparison of total band absorptanoe 
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Table 5.19 Comparison of total band absorptance for the 
2.7)i H 2 O band at T=833 K 
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5.2.4 1.87m B 2 0 Band 


The spectral transmittance results, as calculated by the QRB model, 
are compared with experimental results in Figs. 5.24 and 5.25 for two 
different path lengths but for the same pressure of P e = 765 mm Hg. 

The agreement between these results is quite good for the entire spec- 
tral range of the band except In the wing regions. 

The total band absorptance results calculated by the QRB model 
(only at 300 K) and wide-band correlations are compared with the experi- 
mental results for four different temperatures in Tables 5.21-5.24. It 
is noted that the QRB model can be used to evaluate the band absorptance 
with reasonable accuracy. However „ the agreement between the experi- 
mental results and results of three correlations is seen to be poor, 
especially at low temperatures. 

5.2.5 Total Emissivity for H 2 0 


The total emissivity of water vapor, as calculated by the QRB model 
(only at 300 K) and wide-band correlations, are compared with the ex- 
perimental results of Hottel et al. [9] in Figs. 5.26, 5.28, 5.30 and 
5.32 for four different temperatures and P =* 1 atm. The QRB model re- 
sults do not agree very well with the experimental results at path 
lengths less than 2 atm-cm. However, there is a reasonable agreement 
between these results above 2 atm-cm. Since rotational band contributes 
about 80 percent of the total emissivity at room temperature, the dis- 
crepancy may be due to an error in the description of the line structure 





Table 5.21 Comparison of total band absorptance for the 
1.87|i H2O band at T=300 K 












IP able 5.23 Comparison of total band absorptance for the 
1.87p H2O band at T=833 K 
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in the rotational band. It is noted that the Felske and Tien's results 
have a reasonable agreement with the experimental results for small path 
lengths; however, the Edwards and Balakrishnan, and Tien and Lowder re- 
sults agree well with the experimental results for large path lengths 
except for the high temperature of 2000 K. It should be pointed out 
that although the agreement between the Felske and Tien's results and 
the experimental results is excellent at 2000 K, it does not mean that 
this correlation can estimate the total emissivity of water vapor accu- 
rately at high temperatures. This is because the data from the Hottel’s 
charts, at high temperatures, are estimated from extrapolation and, 
therefore, are less accurate for the entire path lengths. 

Figures 5.27, 5.29, 5.31 and 5.33 show the comparisons of the re- 
sults of total emissivity, as obtained by wide-band correlations, for 
four different temperatures and five different pressures. It is again 
found that the Felske and Tien's results are lower than the results of 
the other two correlations for all cases; and the Tien and Lowder ' s re- 
sults do not have the trend of other results at low pressure. The 
Edwards and Balakrishnan's correlation yields lower values at low pres- 
sures and higher values at high pressures than the Tien and Lowder' s 
correlation. 


5.3 Homogeneous Results for C0 2 + ^0 
Table 5.25 shows the comparison of numerical results of total emis- 
sivity, as calculated by the wide-band correlations for the mixtures of 
water vapor, carbon dioxide and nitrogen, with the experimental results 





TOTAL EKCSSIVITY 








H 2 0 ffOTAL emissivity 





el 







h 2 o total emissivity 












H2O TOTAL EMISSIVITY 




HpO TOTAL EMISSIVTTY 




Table 5.25 Comparison of total gas emiseivities for mixtures 
of water vapor, carbon dioxide and nitrogen 





of Ref. 53. Values of emissivities are corrected for partial overlap- 
ping in accordance with Renner and Varanasi [32]. The Edwards and 
Balakrishnan, and Tien and Lowder correlations overestimate the values, 
while the Felske and Tien's correlation underestimates them. However, 
all these results agree with the experimental results reasonably well. 

5.4 Nonhomogeneous Results 

Since no experimental results of nonhomogeneous gas emissivity are 
available in the literature, comparisons of nonisothermal band absorp- 
tance and emission are made only in this section. 

The results of nonisothermal band absorptance for three different 

bands of CO 2 , as calculated from Eqs. (2.9) to (2.11) by employing the 

wide-band correlations, are compared with the experimental results of 

Edwards et al. [19] in Table 5.26. It is noted that the Edwards and 

Balakrishnan and Tien arid Lowder results have a reasonable agreement 

with the measured values. The Felske and Tien's results are generally 

lower than the experimental data. 

♦ 

Table 5.27 shows the comparisons between the nunerical results of 
nonisothermal band absorptance and experimental results of Edwards et 
al. [19, 54] for two 1^0 bands (6.3 y and 2.7 y ). The results of all 
three correlations show reasonable agreement with the experimental 
results. 

Table 5.28 shows the comparisons between the nonisothermal band 
emission results of 2.7y H 2 O, as calculated with Eq. (2.12) by using 
the wide-band correlations, and the experimental results of Edwards et 
al. [54] and Simmons [55]. The results of correlations, in general, 
show a good agreement with the experimental results except the Felske 











Table 5*27 Comparison of nonisotbermal band absorptance 
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for 2. 7n HoO band 






n i 

H- 

m 

VO 

VO 

CO 

cn 

• 




PM 

? 

'I 

m 

1 

CM 

1 

7 

CM 

1 






c- 






■ 

<U 




VO 

MO 

co 









<T\ 

E'- 

CM 

VO 

m 

vo 

H 

03 

(4 


M 

rH 

en 

-*• 

in 

VO 

c n 

V 

H eg iD 


O 

0 

O 

rH 

rH 

H 

Pi 

a) 

■H 



• 

• 

• 

a 

a 

a j 

+» 

01 

ft 

GH 



O 

0 

O 

O 

O 

O 

r V 




• 







0 




p 

H" 

H“ 

in 


rH 

H 

0) ' 




PM 

rt 

CM 

rl 

O 


. 

0] 




1 

+ 

1 


1 

1 

CM 

t 











+* 













Pi 



’Sf* 

t- 

CO 






<D 



in 

m 

VO 

O 

CM 

H* 


£ 



W 

CM 

in 

in 

rH 

CM 

in 

+» 

a> 



O 

O 

O' 

CM 

CM 

CM 

. W 

•H<HO 



• 

• 

• 

a 

a 

a 


6H 

hi 



O 

0 

0 

O 

O 

O 







■ 





«k 

03 




P 

«sj- 

in • 

H- 

m 

H- 

rH 

• +* 


g 

• 

Pi 

H 

<n 

rH 





’ 

0 

p 


1 

+ 

1 

+ 

+ 

+ 

03 


,a 











03 









« 

01 

*8 

71 

• 


««*• 

m 






■a 



m . 

GO 

t— 

co 

tn 

00 

0 

a) 


n 

CM’ 

in 

in 

rH 

<n 

in 

•H 

01 

•go-a 


O 

• 

O 

• 

O 

• 

CM 

a 

CM 

• 

CM 

a 

03 

•H 

w 

p? 



O 

O 

O 

O 

. 0 

O 

a 

M 





in 



■ 




• 




t- 

CM 

in 




s 

cm 



ej 

cn 

cn 

VO 

H 

in 

VO 

d . 

K 



w 

CM 

h* 

VO 

H 

CM 

in 

.« • 

W 



O 

0 

0 

CM 

CM 

CM 





, 

• 

a 

0 

a 

• 

a 






O 

O 

0 

O 

0 

O 

• 

«H 




» 

H* 



in. 








in 



in 



* « 

0 » 

■as 






• 

CO 

t- 

VO 






. t> 

in 

cn 

CM 

H 

CO 

■H 03 




a 


H - 

VO 

tn 

in 

c- 

+» 03 
H <U 



Pi+» 

* 

a 

a 

a 

• 

• 




m 

0 

O 

O 

O 

O 

0 

td H 
PM Pi 




w 

• m 


co- 









i-i 

H 

rn 

H* 

co 

in 






t 

1 

I 

1 

1 

1 

G 





cr\ 

CM 

CM 

CM 

CM 

in 

3 






H 

H 

rH 

CM 

CM 








t 

rH 

H 

rH 











rH 


PD=((E-E E )/E E )X100 


104 


and Tien's results which are consistently lower than the measured data. 

The discrepancies between numerical and experimental results are 
likely due to improperly approximating the temperature profiles. These 
may be also caused by the experimental error and the Inaccuracy of the 
correlation parameters of the isothermal band absorptance. 


Chapter 6 


CONCLUSIONS 

The results were obtained for spectral transmittance, spectral 
emlssivity, band absorptance, and band emissivity for important infrared 
bands of carbon dioxide and water vapor by employing various absorption 
models in order to establish the validity and usefulness of a simplified 
model for atmospheric and engineering applications. For this purpose, 
selected results were obtained by using the 1 ine-by-1 ivie and quasi- 
random band models. Extensive results were obtained by using the wide- 
band model correlations of Edwards and Balakrishnan, Tien and Lowder, 
and Felske and Tien for a range of temperature, pressure, and path 
length conditions. The results are compared with available experimental 
results. The results were also obtained for total gas emissivities to 
further investigate the usefulness of a model under homogeneous and non- 
homogeneous conditions. 

From a comparison of the spectral transmittance, band absorptance 
and total emissivity at room temperature, it is evident that the quasi - 
random band model results are in good agreement with the experimental 
results in most cases under homogeneous conditions. The quasi-random 
band model can also provide accurate results at higher temperatures if 
proper adjustments are made for overlapping of the lines in a band. The 
model provides accurate results also for nonhomogeneous conditions. 

Thus, use of the quasi-random band model is highly recommended for most 
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atmospheric and engineering applications. Unfortunately, the computa- 
tional costs associated with this model can be prohibitive for many 
applications. 

For engineering applications, use of the wide-band model correla- 
tions is highly recommended. These correlations require significantly 
less computational time and, if used properly, yield results with 
reasonable accuracy. The results of this study reveal that use of the 
Edwards and Balakrishnan and Tien and Lowder correlations is justified 
for carbon dioxide and water vapor at moderate and high pressures and 
relatively large path lengths. The Felske and Tien's correlation is 
useful for low pressures and small path lengths. For nonhomogeneous 
applications, use of the Curtis-Godson approximation, in conjunction 
with an appropriate wide-band correlation, is recommended. The Edwards 
and Balakrishnan and/or Tien and Lowder correlations are recommended for 
most engineering applications involving nonhomogeneous conditions. 

For further study, it is suggested to investigate the feasibility 
of the wideband model correlations for other important molecular 
species by comparing the band absorptance and total emissivity results 
with results of the quasi-random band formulation. The quasi-random 
band formulation can be easily extended to a mixture of molecular 
species. Thus, it would be desirable to investigate the applicability 
of a correlation for a gaseous mixture. This information could be very 
useful for applications involving combustion and flames. 
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A-l 

AL 

AVSI 

BEMI 

BIG 

BIGI 

DEL 

EL 

EM 

FR 

FRB 

FRC 

FRL 

FRU 

IR 

JD 

KR 

LE 

LIB 

LIE 


EXPLANATION OF SYMBOLS USED IN COMPUTER PROGRAMS 


Symbols used in program QRB.FOR 

average half -width of the lines of the molecule 

average of the intensities of the lines within one intensity 
decade 

band emissivity 

upper intensity limit for individual intensity decades 
intensity of the strongest line within one interval 
width of an interval 

energy of the lower states for the individual line 

emissivity from lines within one interval 

frequencies of the individual lines 

frequencies at the interval boundaries 

frequencies at the interval centers 

lower frequency boundary of the band 

^ upper frequency boundary of the band 

index which lets the transmittance be computed for the 
unshifted and shifted meshes 

number of adjacent intervals on both sides of an interval 
from which the contribution is taken into account 

number of interval in the band 

total nunber of lines in the spectrum 

serial number of the first line in an interval 

serial number of the last line in an interval 







■j 
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NSI 

PCK 

PL.PL1 

PNTP 

PREC 

RP 

SI 

SSI 

T1.T2 

TEMC 

TEMR 

TRA 

TRD 

VP 

XI, X2 

A -2 

A,B 

ABSO 

BC 

BIN 

BLP 

BRO 


nunber of lines within one decade of an interval 
Planck's function 

pressure path length of the absorber 
pressure referring to NTP (760 mm Hg) 
effective pressure of the absorber 

exponent to account for the temperature-dependence of the 
rotational partition function 

intensities of the individual lines 

sum of the intensities of all the lines within one decade of 
an interval 

weighting factors used in the numerical integration of 
direct and wing contribution,, respectively 

temperature of the absorber 

reference temperature for the molecular spectral parameters 

array of transmittance values containing results from un- 
shifted- and shifted-mesh computations 

contribution to the total transmittance from lines within an 
interval 

factor accounting for the temperature-dependence of the 
vibrational partition function 

abscissa values for the nunerical integration of the direct 
and wing contribution, respectively 

Symbols used in program NONHQM.FQR 
length 

band absorptance 
band centers 

integrated band intensity at TEMC 
line structure parameter at PREC 
self-broadening parameter 
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BWP 

DENS 

HBIN 

HBLP 

HBWP 

PI.CK 

PML 

PRE 

PREC 

PTOT 

TEMC 

SUM 

A-3 

BETA 

F(I) 

G(I) 

OMG 

RTH 

U(I) 


band width parameter at TEMC 
density of the absorber 
scaled integrated intensity 
scaled line width to spacing ratio 
scaled band width parameter 
Planck's function 
mass path length 
partial pressure of the absorber 
effective broadening pressure 
total pressure 
tanperature of the absorber 
band emission 

Symbols used in subroutine WATER. FOR 
line structure parameter at P e = 1 atm 
line-width-to-spacing temperature-variation parameter 
band-intensity temperature-variation parameter 
-band width parameter at T = 100 K 
integrated band intensity at J = 100 K 
the quantity hcv/kT 
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APPEMOIX B 

QUASI-RANDOM BAND MODEL COMPUTER PROGRAM 


C PROGRAM TO CALCULATE BANO EMISSIVITY WITH QUASI-RANDOM BAND MOOEL 
INTEGER X*W 

DIMENSION FRf 10150) *51(101509 *EL (10150 )*FRB( 85) *FRC( 85) *L IB (65 ) 
4*LIE(65)*MP(85)*BIGI(85)*B1G(6*85)*SSI(5*85) »N5I ( 5*05 ) • 

4AVS1 (5*85) vTRG (83) *X 1(26) «Tl(26) *X2(21)*T2(21)*PCK(85)t 
♦PHI9I *TRA(9*170) 

OATA FH/O. I, 0.4, 1. 0,4. 0*10.0, 40. 0,100.0 •400.0*1000.0/ 

DATA FRL*FRU»DEL/340Q. 0*3800*0* 5.0/ 

OATA PNTP*PREC*TNTPtTEHR/ls 0*0. 1*273. 0*296.0/ 

DATA LE*JO/10104*10/ 

OATA RP*AL/ 1.0 *0.07/ 

OATA Xl/0.0*0.001*0. 0015*0.002*0.003, 0.004, 0.005*0.006,0. 008, 

1 0.01*0.015*0.02*0.03*0.04*0.05*0.06*0.08*0. 10*0.15* 

2 0.2 *0.3*0 .4*0. 5*0.6 *0.8* 1.0/ 

DATA Tl/O .0006 *0.0006 *0.0006 *0.0007 *0.001 *0.001 *0.001 *0.0015* 

1 0.002, 0.003, 0.005*0. 008, 0.01. 0.01*0. 01, 0.015*0.02 *0.03, 

2 0.05 *0.08 *0.1 *0.1*0 si *0.15 *0.2 *0.1/ 

DATA X2/-i«0*-0.9»-0.8»— Q.7»-0.6*— 0.5*— 0*4*— 0.3»— 0.2»— 0.1* 

1 0.0*0.10*0.20*0.30,0,40*0.50*0.60*0. 70*0. 80*0.90* 1.0/ 

OATA T2/0. 1 *0 .4 *0.2*0.4* 0.2*0. 4*0.2 *0. 4*0.2 *0.4* 

1 0.2 *0.4*0. 2 *0.4*0. 2 *0.4 *0.2 *0.4 *0.2 *0.4 *0.1/ 

ACCEPT *. TEMC 
ACCEPT 4* VP 

READ(23*13) (FR(X)*SI (X)*ELtX)*X-l*LE) 

13 F0RMAT(F10.4,E15.S,F12.4) ' 

C CONVERTS THE LINE INTENSITIES AND THE AVERAGED LINE-WIDTHS FROM 
C REFERENCE TD AMBIENT CONDITIONS OF TEMPERATURE AND PRESSURE 
CST— ( S.0RTI TENR/TEMC 1 )*PREC/PNTP 
PA RT- VP * ( ( TENR/TEMC |44Rp| *2 .69E+19 
FACT-1.43 9* (TEMC-TEMR )/( TEMC*TEHR I 
DO 101 X-1*LE 

101 SIIX)«SI(XJ4PART4EXP(EL<X)4FACT) 

ALA-AL*CST 

C INITIALIZES THE FREQUENCY INTERVAL BOUNDARIES FOR THE UNSHIFTED MESH 
DELA-0.5*DEL 
RK— (FRU— FRL1/0EL+0.1 
KR-RK 

FR8(1)-FRL 
DO 201 K-1,KR 
FRB(K+1)-FRB(KB*DEL 
FRC(K)-FR8(K)«0ELA 

C CALCULATES PLANCK FUNCTION FOR EACH FREQUENCY INTERVAL 
RNUN-1.192540EL«FRC(K)4FRC(K)4FRC(K)41.E-05 
EEX-1.43S8*FRC(K) 

201 PCK(K)-RNUM/(EXP(EEX/TEMC)-1.) 

WRITE(3*222> PREC«TEMC 

222 FORMAT ( ////10X*16HTHE PRESSURE IS *F10.4*4H ATM* 
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*10X*19HTHE TEMPERATURE IS *F10.2*2H KJ 
IR-0 

C STARTS THE COMPUTATION FQF THE UNSHIFTED MESH 
CO TO 121 
122 KR-KR-1 

FRBUI-FRL+OELA 

C DEFINES THE FREQUENCIES AT THE INTERVAL BOUNDARIES AND CENTERS 
121 DO 100 K-l.KR 

FRB(K*1)-FRHK)*0EL 
100 FRCIK)-FRB(K>*DELA 

C DETERMINES THE NUMBER OF LINES IN EACH FREQUENCY INTERVAL 
X-0 

DO 102 K-l.KR 
H-0 

104 IFf XaGEaLE) CO TO 103 
X-X+l 

IF(FR(XI«LT«FRBIK1) CO TO 104 
IF(FR(XI.GE.FRBtK*lM CO TO 105 
H-M*l 
CO TO 104 

105 IF(K.CE.KR) CO TO 103 
X-X-l 

103 HPIK1-N 
102 CONTINUE 

C ASSIGNS NUMBERS TO THE FIRST AND LAST LINE OF EACH INTERVAL 
X-0 

00 106 K«1*KR 
IFIHP(K).EQ.O) CO TO 106 
LIB(K)-X*1 

LIE(K)-LIB(K)*HP(K)-1 
X-LIE(K) 

106 CONTINUE 

FIND OUT THE HIGHEST INTENSITY IN THE INTERVAL* GENERATES THE 
INTENSITY OECADES*OISTRIBUTES LINES INTO INTENSITY DECADES* 
DETERMINES 1 THE NUMBER OF LINES IN EACH OECAOE AND THEN COMPUTES 
AVERAGE INTENSITY FOR EACH OECADE 
DO 107 K— 1*KR 
IFIMP(K).EQ.O) GO TO 440 
4B-LIBIK1 
JE-LIEfK) 

BIGKKI-SIf JBI 
JC— JB*1 

00 108 J-JC*JE 

IFIBIGI (KI«GE«SI(J) I GO TO 108 
BICKKI-SH J> 

108 CONTINUE 
DO 109 1-1*6 
IX— 1*1 

109 BIC(IoKI-BlGI(Kl*10.**IX 
00 110 1-1*5 
N-0 

SSIf I*K)-0« 

DO 111 J-JB.JE 

IF(SI(J).GT.BIG(I*K1) GO TO 111 
IF(SI(JI.LE.BIC(I*1*K)> CO TO 111 
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N-N*l 

SSI(I,K)-SSICI.KI+SI< J) 

U1 CONTINUE 
NSI tl vKI-N 

IFINSKIvKKCE.il CO TO 1X0 
NSIfI«K)-l 

110 AVSIf I vK)— SSIf I vK)/FLOATf NSI tlvK) I 
CO TO 10? 

440 00 Ml I-lv5 
NSItlvKl-l 

441 AVSIfI*K)-0. 

107 CONTINUE 

PI-3.14159 

RHO-ALA/DELA 

C STARTS THE CONFUTATION OF TRANSNITTANCES FOR EACH INTERVAL 
DO 301 IPL— 1»9 
PL-PUCIPLI 
00 112 K-lvKR 
TRCtKI— 1. 

DO 113 J-lvKR 
TRO-1. 

JA-IABSU-K) 

IFt <IA.CT.JOI CO TO 113 
ZI-FRCfKI—FRCf 4) 

EPS1— ZI/OELA 
00 115 I-lv5 
NSJ-NSKIvJI 

PNUH-RHO*RHO*AVS I C I . J I *PL/ * PI+ALA I 
RES-O. 

IFt J.NE.KI CO TO 116 

C EVALUATES THE INTECRAL FOR DIRECT CONTRIBUTION 
DO 117 M-lv26 

YY-PNUH/tXlU)*XltUI*RHO*RHO> 

IFfYY.CT.90.) CO TO 119 
Y-EXPt-YYI* 

CO TO ,117 

119 Y— 0. 

117 RES-RE$+Y*Tlf II) 

CO TO 244 

C EVALUATES THE INTECRAL FOR MXNC CONTRIBUTION 
116 00 118 W-lv21 

Y Y- PNUH/ < f E PS I~X2 ( H ) ) * < EP S I-X2 1 M 1 1 ) 

IFtYY.GT.90.) CO TO 120 

V-EXPt-YYl 

CO TO 118 

120 V-0. 

118 RES-RES*Y*T2fU) 

RE5-RE5/6. 

244 IFt RES.EQ.O.) CO TO 234 
255 TRD-TR0«RES**NS4 
CO TO 115 
234 TRD-O. 

115 CONTINUE 

113 TRCfK)-TRC(K)*TRO 

112 CONTINUE 
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IFiIR.CE.il CO TO 123 
00 114 K»1«KR 
114 TRA<IFL*2*K-lt«TRCtK) 

CO TO 301 

123 00 124 K-1*KR 

124 TRA(IFL*24K)-TRC(K) 

301 CONTINUE 
IR-IR»1 

C COES TO 122 ANO STARTS CONFUTATION FOR THE SHIFTED NESH 
IFCIR.LT. 21 CO TO 122 

C AVERAGES RESULTS OVER SHIFTEO ANO UNSHIFTED NESHES 
OEN-5.660*ITEHC**4 I41.E—05 
DO 302 I«l .9 
00 125 K-l.KR 

125 TRC f K)»i TRA M » 2*K-1 ) *TRA CI»2*K> *TRA i I *2*K*l 1 1 /3. 
URITE(3*601 (TRCIK)*K-1«KR) 

60 FORMAT (///(10F12.5/I1 
C CALCULATES ENIS5IVITY 
SUM-0.0 
DO 220 K-i*KR 
EH-1.-TRGIK1 

220 SUM«SUM*EM*FCK IK) 

BEMI-FI*SUN/DEN 

URITEf 3*221) BEHI.FLltll 

221 FORMAT ( // 10X* 10HTHE EMISSIVITV IS *E15.7* 

* lOXtlSHTHE FATHLENCTH IS *F10.2«?H ATH-CMI 

302 CONTINUE 
STOF 
END 
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APPENDIX C 

COMPUTER PROGRAM TO CALCULATE NONHOMOGENEOUS 
BAND ABSORPTION ANO BAND EMISSION 


C MAIN PROGRAM TO CALCULATE NONHOMOGENEOUS TOTAL BAND ABSORPTANCE ANO 
C BAND EMISSION 

DIMENSION PLCK(6*290l*ABSO(6*29OI*BIN(8) *BLP (8) «BUP( B I *HBIN(6) • 

* HBLPI6>*HBWP<6)*HI6IN<6>*HIBLF(6>»HIBMP(6)*SUN(6>*BCt6) 

TlfY)-3.92*V«1019.0 
C T2IYI— 7.33*V*1439.B 

DATA BC/lAO.a*9O0.O«ia0O.O«3760.O«53SO < ,O v 7250.O/ 

A-0.0 
8-60.0 
C C-60.0 

XN1-120.0 
C XN2-60.0 

X1-(B-A>/XN1 
C X2-(C-B)/XN2 

PTOT-I.O 
PRE-0.786 
00 90 1 - 1*6 
MBINIIJ-0.0 
HBLPIH-0.0 
90 HBUP(I)-0.0 
PML-0.0 
M-l 

400'ARG-A+FLOAT(MI*X1-X1/2. 

IFIARG.GT.BI GO TO 200 
ARG-BMFL0ATIMI«»XNll*X2— X2/2. 

IF(ARG.GT.t) GO TO 200 
0LTY-X2 
TEMC-T2IARG} 

TEMC2-T2I ARG-X2/2. I 
GO TO 300 
100 QLTY-X1 

TENC-TKARG) 

TEMC2-T1I ARG—X1/2. ) 

C 300 PRE-PRESSURG) 

300 DENS-2193. 119+PRE/TEMC 

CALL HATER (TEMC*BIN*BLP* BMP) 

PML-PML*DENS*DLTY 
BR0-8.6+SQRT ( 100.0/TEMC) +0.9 
P REC-PTOT +BRO*PRE 
DO 150 K-l*6 

HBIN(KI-HBIN(K)«BINfK)*DENS*DLTY 
HBLP(K)-HBLP(Kl*BLP(K)*PREC*filN(K>*OENS*DLTY 
HBMP f K J— HBHP ( K ) +BWP ( K ) *B I N( K ) *0ENS*0LTY 
HIBINIKI-HBINUI/PML 
HlBLP(K)-HBLP(K)/(PML*HIBtNf Kl ) 

HIBNP(KI-HBUPfK)/(PML«HIBIN(Kl) 


original page ra 
OF POOR QUALITY, 


A8SO(K»N)-EOUAROCPNI.tl.*HlBINm*HI8LP(K)«HIflUP(Kt> 
PUX(K«NI-PLANCK(TEHC2»BC(KD 
ISO CONTINUE 

MRlTEOtlOl CABSOCI. Ml. 1-1,61 
10 F0RHATI/6F19.7I 
H-H*l 
CO TO 400 
200 00 290 J-l»6 

HBINt JI-HBZN( JI/PItL 
HBLPI JI»HBLPU1/(PHL*HBIN(4I I 
HBWPI JI-HBMPC JI/IPHL*HBINU)I 
ABS0RP-E0UAR0(PNLtl**H8IN(JI»HBLP(J).HBUP(J)> 
WRITE(3,2Q) PNL,HBIN(J),HBLP( J),HBMP( J) tABSQRP 
20 FORftATf // *9F16«7> 

290 CONTINUE 

00 900 N-l,6 
900 SUIHNI-O.O 

00 600 L-ltN-2 
00 700 1-1,6 

700 SUN (I)-SUH( 1 I-HA8S0C I ,L+1 >-ABSO( i.L )) *PLCK( K *L*1 J 
MRITE(3,30| (;UN(II *1-1*6) 

30 F0RRAT<76E19.7I 
600 CONTINUE 
STOP 
END 
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SUBROUTINE WATER (TEMC*BIN*BLP*BUP1 

DIMENSION 8ETA(81*RTH(8l*QMGf S1*BIN(8I*8LP(81*BWP(61«U(31 « 

* SUMI (3 1 *HEXf 31 * EXHOl *F (61 
C BAND ABSORPTION PARAMETERS 

DATA RTH/44209. 0*9200. 0*41. ZtO. 19*3. O f 2. 9*2. 3o22. 4/ 

DATA BETA/O. 14311*0. 1431 1*0 .09427*0. 13219*0. 0B169t0. 1162 8* 

* 0.13219*0.13219/ 

DATA 0M6/69 .3*28. 4*96.4*60.0*43. 1*32. 0*60. 0*60.0/ 

C CALCULATE TEMPERATURE OEPENDENTvFUNCTION FOR LINE WIDTH PARAMETER 
U(ll-9294.9/TEMC 
UC2 1-2294. 9/TEMC 
U(3>— 9404.1/7EMC 
DO 100 1-1*3 
SUM-0.0 
V-0.0 

1 ADD-SQRT(¥*l. I*EXP (-U( Il*¥/2.1 
IF! ADD.LT.1.E-04I 60 TO 2 
SUH-3UH+ADD 

V-V+l. 

60 TO 1 

2 SUM I( I l-SUM*SUM 
100 CONTINUE 

SUH-0.0 
¥«fc. 0 

3 AOD-SORTI V+l. 14SQRT(¥-»2. )*EXPl-U(2l«V/2. ) 

IF(A00.LT.l.E-04) 60 TO 4 

SUM- SUH+ADD 
¥-¥♦1. 

60 TO 3 

4 SUH1I-SUM+SUH 
DO 200 1-1*3 
EXN(I1-1.-EXP(-U(I11 

200 HEX(n-ll.-EXP(-Ufn/2.1)*<l.-EXP(-U(l 1/2.1 1 
C 

Ffll-EXMfri*EXM(2)+EXMI2)*EXH(3l*SUHX(2>/<HEX<l)*HEX<3)> 

F(2 l-EXMf 1I*(EXMC 2l**3l*EXM(3>*SUMII/f 2.*HEXC1I*HEX<3II 
F(3!-E*XM( 114EXH(2>4EXM(ZI*EXH(31*EXH(31*SUMI J21*5UNI(31/HEX(11 
FI4)-EXM(11 *EXH(11*EXM(21*EXM(31*EXN(3 1*SUNI (11*SUMI < 31 /HEX (2) 
F(31-EXH(11*EXM(11*EXM(21*EXH(31*$UHI (11/(HEX(2)*HEX(31 1 
Ff6l-EXM«ll4EXMC2l*EXM(3»*EXM(3f*SUMIC3l/(HEXCl>*HEXC2)) 

C 

TEM-SORTITEHC/IOO.OI 
C CALCULATE BAND PARAMETER 

B IN (1) -RTH 1 1 1 *EXP I -9 . 0/T EM 1 
BWP(1I-QHG(1)*TEN 
BLP <11 -BETA 111 /TEN 
BIN <2 1 -RTH 1 2 )4£XP<— 17.6/TEM1 
BWP(21— OHG( 21*TEM 
BLP (21-BETA (2 ) /TEM 
C 

RTH (4 1-0. 194 ( l.-EXP l-2.*U (2 1 1 1/ (EXH ( 2 1 4£XM(2 ) 1 
RTH 1 9 1-3. 04 1 1. -EXP I— I U (2 1+U (31 1 1 1/(EXM(214EXM(31 ) 

• RTH (61-2.94 (l.-EXP (-(U(l 1+U (3) 1 1 1/(EXN(114£V.h(3) ) 

DO 300 1-1*6 
BIN (I +21— RTH (I +21 
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BMMI*2)»0nG(I«2)*rEN 
300 BLP (I*2)"BETA(I*2) +F( I J/TEM 
<4, * elN,7,+BINC 

?. N »7"? ORT(RTHM,#B L» , Mn*SQRTIRTH«7l 

BLP<4I«XNUH*XNUH/B1N(41 

RETURN 


♦ BLP 171 >+5QRT(RTH(8) *BLP 
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SUBROUTINE CARBON! TENC.B IN*BLP,BHP) 

DIMENSION BETA (6 ) «RTH(6) , OMG( 6) ,BIN(6 ) ,BLP(6 1 »BUP (6 ) *U(3 ) tSUMI 13)* 
■> SUNI I (21 ,EXM(3) *HEX(3) *F(6)*G(6)«H(6) 

DATA RTH/19.0,2.47E-09*2.48E-09*110. 0*4. 0*0.66/ 

OATA BETA/0. 06157,0. 04017, 0.11886,0. 24723*0. 13341,0.39305/ 

DATA OMC/12.7,13.4,10.1,11.2,23.5,34.5/ 

UU1-1943.82/TEMC 

U(2)-939.66/TEHC 

U ( 3 I— 3447. 365/TEMC 

DO 100 1-1*3 

SUM-0.0 

V-0.0 

1 AOO-SQRT ( V*l. )*EXP<-U( I)*V/2. ) 

IFI ADD.LE.l.E-041 CO TO 2 
SUM-SUM 4-A00 . 

V-V*l. 

CO TO 1 

2 SUMI(I)-SUH+SUH 
100 CONTINUE 

DO 200 1-1*2 

SUM-0.0 

V-0.0 

3 AD0-5QRT(V+l.)*SQRT(V+2. l*EXPf-U(I»*V/2.» 

IF( ADD.LE.l.E-04) GO TO 4 
SUM-SUH+ADD 

V-V+l. 

GO TO 3 

4 SUMIK 1 )— SUH*SUH 
200 CONTINUE 

SUI-0.0 

V-1.0 

5 ADD-SORTC V+1.I*EXP(— UHM'V/2. ) 

IF( ADD.LE.l.E-04) CO TO 6 
SUI-SUI+ADO 

V-V+i. 

GO TO $ 

6 SUI-SUI+SUI 
00 300 1-1*3 

EXM (I)-l.-EXP <-U(I > > 

300 f\'EXIII»Cl.-EXP(-U(II/2.)»6Cl # -EXP(-Um/2.l) 

UENF— W(EXM(l)*EXM(in-l. 

c 

F(l)— EXH( 2) *EXM(2I*EXN< 1 I *EXH( 2)*EXM( 3I*SUNI I (2) / (2. *HEX ( D4HEX (3 ) 

♦i 

FC2I-EXM(2I*EXM(2>*EXM(3I+6XM(3)+EXM(3»*SUI+SUMI(2)*SUMI(3)/DENF 

F(3)-F(2) 

F(4l— EXH( 1) *EXN(2)*EXN(2 I *EXH(3I*EXH (3 )*SUMI (2)*SUMI ( 3)/HEX( 1) 
F(5)-EXH(1) *EXN(1)*EXM(2 I 4EXM( 2I*EXM(3 >*EXM( 31+SUMI (11+SUNK2I+SUM 
41(3) 

F(6)-EXM(1)*EXR(1)*EXM(1)*EXM(2)«EXM(2)*EXM(3)4EXM(3)4SUMII(1)«SUM 
*1(2 l*SUHI (3 )/2. 

C 

G(l)> 1.0 

C(2 )“1.383*(1.-EXP(U(1I— U(3I I I* (2 .—EXP (**U(1) ) l*EXP( — U (1) )/(EXM(l)4 
4EXMI3 ) )*1.EQ8 


f 
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G(3>-1.303*U.-EXP<2.*U(2>-U(311)*(2.-EXP(-Um>J*EXP (-U( 11 >/(EXH< 
*l)*EKn(3>>*l.E03 
6(41-1.0 

e(5)-(i.**EXP(-u(i}-uocn/iExnm«Exnm» 

C(6I-a.-EXP(-2.*um-U(3m/(EXn(ll4EXN(l)*EXn(3ll 

Hill— 1.0291 
H(2)-1.0273 
H (3 >-1.0273 
H(4>— 1.0237 
H(5>— 1.0238 
H(6>— 1.0237 

TEB-SQRTCTEKC/IOO.OI 
00 400 R-l*6 
BINOl)-RTH(N)*6(N> 

tolp('a)-b«t«(H)*r (■J/(t««*h(a)> 

400 BWPf HI— OHG(HJ+TEH 
RETURN 
END 


